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A modified version of the spinless Anderson model is studied by means of the continuous¬ 
time quantum Monte Carlo method. This study is motivated by the peculiar heavy-fermion 
behavior observed in certain Samarium compounds, which is insensitive to magnetic field. The 
model involves M channels for conduction electrons, all of which interact with local / electron 
via the Coulomb repulsion Ufc, while only one channel has hybridization with the local state. 

The effective hybridization is reduced by the Anderson orthogonality effect, and a quantum 
critical point occurs with increasing M and/or increasing Ufa- The numerical results at finite 
temperature of the local charge susceptibility are well fitted by a simple scaling theory for 
all M. However, the single-particle spectrum is described by a double Lorentzian for M > 1, 
in contrast with the single Lorentzian with M = 1. A quasi-particle perturbation theory is 
presented that reproduces the quantum critical point for large M. The quasi-particle theory 
gives not only the renormalized energy scale, but its extrapolation toward higher energies being 
consistent with the double Lorentzian spectrum. 

KEYWORDS: charge Kondo effect, continuous-time quantum Monte Carlo method, quasi-particle pertur¬ 
bation theory, thermodynamic and dynamic properties 


1. Introduction 

Recently, peculiar heavy-fermion behavior has at¬ 
tracted attention in certain Samarium compounds with 
large specific heat coefficient 7 which is insensitive to 
external magnetic field. For example, the filled skutteru- 
dite compound SmOs4Sbi2 has 7 ~ 0 . 8 J/(K^ • mol) even 
though it is mixed valent.Similar behavior has been 
found in systems such as SmPt4Gei2^^ and SmT2Al2o 
with T=Ti,V,Cr,Ta.^“®^ The resistivity of SmT2Al2o 
shows clear Kondo-like logarithmic temperature depen¬ 
dence, which however is insensitive to external magnetic 
field.It has been suspected that charge degrees of free¬ 
dom is responsible for the heavy mass because of the 
field-insensitivity, in striking contrast to ordinary Kondo 
effect which is sensitive to magnetic field. 

Motivated by these experimental observations we 
search for a charge fluctuation mechanism that gives 
rise to energy scale much smaller than bare hybridiza¬ 
tion. As the simplest attempt, we study the (spinless) 
multichannel interacting resonant level model (MIRLM) 
by means of the continuous-time quantum Monte Carlo 
method, which starts from the Anderson model involving 
the hybridization {V) of the local charge with only one of 
the conduction electron orbitals, but includes additional 
Coulomb interaction {Ufc) felt by all conduction orbitals 
with the local / state. The Hamiltonian of this model is 
written as 

'H ~ He + Hf -\- Hhyb + Hfc 


M-l 

= EE £kcl fCk,£ + Efpf -f V(pco + cj/) 

k £=0 

(1) 

where M is the number of the conduction electron chan¬ 
nels, Cfc7 (c][, ^) is the annihilation (creation) operator 
of the Bloch state k in the .^th channel, and C£ = 
jy-1/2 with N being the number of sites. 

The model given by Eq. ( 1 ) leads to rich physics under 
finite values of the Coulomb interaction Ufc with more 
than one conduction channels (M > 1 ). This is due to 
the increasing dominance of the Anderson orthogonality 
effect arising from the screening channels over the exci- 
ton enhancement coming from the hybridizing channel. 
Namely, the presence of multiple channels of conduction 
electrons leads to non-trivial low-energy renormalization 
of bare V. 

In addition to the motivation provided by the peculiar 
heavy-fermion state of certain Samarium compounds, 
such quantum impurity models play important role to 
understand many-body phenomena realized in single ar¬ 
tificial atoms, or quantum dots, with multiple levels in¬ 
teracting with a Fermi sea and encountering further local 
interactions. 

The single-channel version of the model has been ex¬ 
tensively studied by many authors with various ana¬ 
lytic methods including bosonization,®) Bethe ansatz,^^ 
Anderson-Yuval mapping to Coulomb gas®^ and pertur- 
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bative renormalization.®’ However, much less is known 
about the multichannel version of the model. Although 
the original interacting resonant level model has been ex¬ 
tended to multiple channels by perturbative renormaliza¬ 
tion^^^ and also been studied by numerical renormaliza¬ 
tion group,these works did not discuss the dynam¬ 
ics under finite hybridization and at finite temperatures. 

In a previous paper,we already studied numeri¬ 
cally the single-channel (M = 1) version of the model 
in the negative Ufc range and provided quantitative in¬ 
formation about the dynamics at finite temperatures. 
In this paper we extend the study for the positive Ufa 
range and for the presence of multiple conduction chan¬ 
nels. By means of the continuous-time quantum Monte 
Carlo method we investigate both thermodynamic and 
dynamic properties of MIRLM in a wide range of the 
Coulomb interaction Ufc- Especially, we are interested 
in the applicability and accuracy of perturbative ap¬ 
proaches. 

This paper is organized as follows. In Section 2 we de¬ 
rive the renormalized hybridization for the multi-channel 
case by means of perturbative renormalization within a 
simple scaling theory. In Section 3 the continuous-time 
quantum Monte Carlo algorithm is formulated. The nu¬ 
merically obtained static and dynamic properties are pre¬ 
sented in Sections 4 and 5. We construct a quasi-particle 
perturbation theory in Section 6 in order to better under¬ 
stand the numerical results. Finally, Section 7 is devoted 
to the summary of this paper. 


2. Perturbative Renormalization Approach 

2.1 Scaling energy with multiple conduction channels 
The multichannel version of the interacting resonant 
level model was first introduced in Ref.ll, which mapped 
the MIRLM to an anisotropic Kondo-like model and de¬ 
rived the scaling equations of the corresponding Kondo 
model. Here, we follow a different way. Namely, we ex¬ 
tend the method that we used in our previous work^^^ 
to the case of multi channels for the conduction elec¬ 
trons to obtain the hybridization renormalized by the 
Coulomb interaction. In that work we derived the effec¬ 
tive hybridization for the single-channel case by using the 
effective Hamiltonian method^®^ to take account of the 
simultaneous effect of the hybridization V and Coulomb 
interaction Ufc- If we include multiple conduction chan¬ 
nels, the expression for the renormalized hybridization is 
modified as 


, , 5D I , 6D 

V' = V{l-u- + -u^M- 


( 2 ) 


where D is the band cutoff and we introduced a dimen¬ 
sionless coupling constant u = poUfc with po being the 
density of states. The factor M accounts for the closed 
conduction-electron loops that represent the Anderson 
orthogonality effect. Following the same procedure as in 
Ref. 14, the effective hybridization is finally obtained as 

V*=vi-^\ , (3) 


where Aq = irpoV^ is introduced. Correspondingly, V* 
defines the scaling energy A*^ = Trpo{V*)‘^, and we ex¬ 


press from Eq. (3) as 

A* A* 

with 

rj = u(2 — Mu). (5) 

Coefficient rj controls A*^, and therefore the effective hy¬ 
bridization V* . Namely, the exciton effect (~ u) com¬ 
ing from the single hybridizing channel enhances the hy¬ 
bridization with increasing Coulomb interaction, while 
the Anderson orthogonality effect (~ —Mu^) from all 
channels acts against the exciton effect by blocking the 
hybridization. As a result, the competition of these two 
effects renormalizes the hybridization at low energies in 
a non-trivial way. 

The renormalized hybridization obtained in Eq. (3) is 
valid only for small values of the bare parameters V and 
Ufc because of the perturbative renormalization treat¬ 
ment. Making use of the analogy with the x-ray threshold 
problem^®) by considering finite Ufc but neglecting its in¬ 
terference effect with the infinitesimal hybridization, an 
associated phase shift can be introduced®’ as 

u = — = —tan“^(7rpoC7/c/2) (6) 

TT TT 

in place oi u = poUfc to account for the multiple scat¬ 
tering by Ufc to infinite order. 

We have found in Ref. 14 that the phase shift picture 
accounts quite accurately for the effective hybridization 
in the negative Ufc range in the single-channel case. 


2.2 Scaling energy at finite temperature 

Let us first quote Schlottmann’s extension^®) of the en¬ 
ergy scale A*g to finite temperatures (T) and frequencies 
(e); 


In^^^ =77 
An 


log 


D 

2^ 



AL(£,T) 

2-kT 



( 7 ) 


where tp is the digamma function. Regarding the energy 
dependence of A*^(e, T), the above expression is an inter¬ 
polation formula between A*^ given in Eq. (4) at £ = 0, 
and the result of the x-ray edge problem Ao(|e|/D)’' for 
|£| ^ A*^. Based on the idea that the interacting reso¬ 
nant level problem is still described by a resonance with 
the renormalized width A*^(£, T) under finite Ufc instead 
of the bare Ag, Schlottmann obtained the charge suscep¬ 
tibility as^®) 


x:{e,T) = 


1 

7 r£ (£ -k i2A*^{T)) 


X 



^tciT) , e \ 
2ttT 2ttT ) 


- V' 


Kc{T) \ 

2-kT ) 


( 8 ) 


by taking the convolution of two simple resonances as¬ 
sociated to the /-electron Green’s function. The static 
component is given as Xo = = 0), which is obtained 
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as 


Xo(T = 0 ) 


1 


(9) 


at zero temperature from Eq. ( 8 ). 

Equation (9) expresses that the charge susceptibility 
is scaled with the single scaling energy A*^. This can be 
understood by recalling the Ward identity,which 
is a constraint for the correlation functions dictated by 
conservation laws present in a model. In the case of the 
resonant level problem the total number of the local and 
conduction electrons, i.e. the charge is conserved, which 
leads to a relation between the vertex corrections and self 
energies.As it was shown by Schlottmann in Ref. 10, 
a consequence of the Ward identity for the resonant level 
problem is that the vertex corrections and self-energy 
cancel each other in quasi-particles, and thus the quasi¬ 
particle density of states pf{e) is completely determined 
by the scaling energy A*^ as 


P/(e) 


al 


7r((e-e/)2 + (A*J2)- 


( 10 ) 


By further use of the Ward identity, the following relation 
was obtained^®^ between the charge susceptibility and 
specific heat coefficient 7 at zero temperature: 

Sy/TT^ = X* = = -^, (11) 


i.e. Xo is entirely determined by non-interacting quasi¬ 
particles through the resonant level given in Eq. (10). 
We can apply the argument above also to the case with 
M > 1. 


2.3 Quantum critical points 

The competition of the exciton effect with the Ander¬ 
son orthogonality effect reflected in coefficient rj given 
in Eq. (5) drives the system toward a quantum critical 
point at about u ^ 2/M. To be more precise, we write 
the exponent x{u, M) in Eq. (3) as 


with 


x(u, M) 


—u + Mu^/2 
1 + 2u — Mu^ 

1 1 
2 2M{u — u-){u — u+) 


u±{M) 


±(i±VTTm). 


( 12 ) 


(13) 


When u approaches to u±, the exponent x{u, M) di¬ 
verges to positive infinity, which means that the effec¬ 
tive hybridization V* vanishes, i.e. the local / charge 
becomes decoupled from the conduction electrons. The 
point R* = 0 is equivalent with the ferromagnetic Kondo 
fixed point with degeneracy between the empty and oc¬ 
cupied / states. We note that the value |u|± ~ 0(1/'/M) 
obtained in Eq. (13) for M ^ 1 is within the range of 
the perturbative treatment, while |it|± ~ 0{1) for M = 1 
might be artificial. 

Taking the positive Ufc range, the perturbative treat¬ 
ment with u = poUfc predicts a quantum critical point 
at It = u+ = 1 -|- even for the single-channel case. On 



Fig. 1. (Color online) An example of the Monte Carlo configura¬ 
tion. 


the other hand, in the phase shift picture, the condition 
1 -|- 2u — = 0 of vanishing hybridization requires 

M > 3 with u = 1 corresponding to the maximum phase 
shift Sjj = IT. The study of MIRLM by means of numer¬ 
ical renormalization group method^^’ found a satura¬ 
tion of the renormalized hybridization for M = 1 in the 
positive Ufc range, and vanishing V* with increasing Ufc 
for M > 3. However, it remains to see what happens for 
M = 2. Although the NRG study found a suppression 
of the renormalized hybridization for M = 2 by increas¬ 
ing Ufc, it cannot be decided explicitly from the data 
whether it reaches zero or not. 

3. Continuous-Time Quantum Monte Carlo Ap¬ 
proach 

In this section, we analyze the MIRLM using the 
continuous-time quantum Monte Carlo method. In 
the previous paper,we presented an algorithm for 
the single-channel case, M = 1, based on an expan¬ 
sion with respect to V and Ufc. The advantage of the 
double-expansion algorithm compared with an ordinary 
weak-coupling expansion with respect to Ufc is that the 
computational cost increases only linearly as M is in¬ 
creased, yielding efficient calculations for large M. Since 
the extension to M > 1 is straightforward, we only briefly 
describe difference from the case with M = 1 in the fol¬ 
lowing. 

The Monte Carlo configuration is expressed schemat¬ 
ically in Eig. 1. Here the imaginary time sequence r = 
{ti, t{, ..., Tq, Tg} defines the R-expansion process. The 
interaction Ufc is described by a time-dependent poten¬ 
tial which fluctuates between uq and ui depending on 
the occupation of the / state.The existence of the 
extra channels M > 1 modifies the weight Wc, which is 
factorized as 

M-l 

Wc(r,po,--- ,pm-i) = detT)o(T,Ro) detidi(p^), (14) 

e=i 

where fii denotes a set of imaginary times at which scat¬ 
tering takes place. Dg is a mi x mi matrix consisting of 
the bare Green’s functions connecting two time points in 
p.1- 

The last factor det requires an extra cost compared 
to the case of M = 1. The important point is that the 
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Fig. 2. (Color online) Histograms Pq and PrriQ for different values 
of channel number M. Here q denotes the ^-expansion order and 
mo the {/j-(.-expansion order in the 0th channel. The parameter 
values are chosen as Ufa = 0.4, V = 0.2 and /3 = 200. 

expansion order mi for each channel is not changed so 
much when M is increased, as confirmed from the his¬ 
togram in Fig. 2. It means that the size of the matrix 
D( is almost independent of M and therefore, the com¬ 
putational cost increases only linearly against M. In the 
weak-coupling algorithm, on the other hand, one com¬ 
putes the determinant of a mi x mi matrix for the 
/ state. It means that the matrix size is proportional to 
M and hence, the cost increases according to 

The single-particle Green’s function consists of three 
components, Gf,Gc and G/c. The relation between each 
component and the self-energy is summarized in Ap¬ 
pendix A. In numerical calculations, we use the con¬ 
stant density of states po = 1/{2D) with the band cut¬ 
off U = 1 for all channels. We restrict ourselves to the 
particle-hole symmetric case, £/ = 0. Spectra in the 
real-frequency domain are obtained by analytic contin¬ 
uation iSn e + iS using the Fade approximation. We 
imposed the condition for the particle-hole symmetry, 
ReG/(f£„) = KeGciiSn) = InrG/c(fen) = 0, to improve 
accuracy. 

4. Static Properties 

We first discuss the static charge susceptibility Xoi'T) 
to check applicability of the scaling theory in Section 2. 
By taking the limit e = 0 in Eq. (8), we obtain the ana¬ 
lytic expression Xo('^) static charge susceptibility 

as 

= + (15) 

In comparing this expression with our numerical results, 
we should take care of influence of the band cutoff D: 
the analytic expression was derived in the limit D ^ oo, 
while the simulation is performed with a finite value of 
D. We introduce a correction factor a = 1 — 2A/{ttD) 
and replace Xg with cxxg to take the influence of finite D 
into account. 

Top part of Fig. 3 shows comparison between the nu¬ 
merical results xo and the analytic expression axg for 
M = 1 and M = 5. The analytic expression turns out 
to give an excellent fit of the numerical data in the wide 




Fig. 3. (Color online) Numerical result for static charge suscepti¬ 
bility as a function of f/j-j. for /3 = 200. Top; Comparison between 
M = 1 and M = 5 with V = 0.2 fixed. Bottom: Comparison 
of different values of V for M = 1. The inset shows the Uf^- 
dependence of the zero-energy limit of the imaginary part of the 
charge susceptibility. Schlottmann’s susceptibility Xq is shown 
by solid green lines. 

range of Ufc- Thus, we conclude that xo is scaled with a 
single parameter A*^, as it was discussed in Section 2. 
We also make a comparison of the imaginary part of 
the charge susceptibility, lmxc(e)/e|e= 0 j in the inset of 
Fig. 3. It can be seen that the agreement is good for 
moderate values of Ufc- However, a systematic deviation 
is observed for a large value of Ufc with M = 1. It might 
be caused by the neglect of the energy dependence of A*^, 
in Schlottmann’s formula (8). 

Finally, we comment on the hybridization-dependence 
of the charge susceptibility. Although arbitrary value of 
the Coulomb interaction Ufc can be handled by the as¬ 
sociated phase shift introduced in Eq. (6), Schlottmann’s 
formula (15) is valid only for small values of the bare hy¬ 
bridization. Thus, it is interesting to check the validity of 
Xo as the value of the bare hybridization V is increased. 
Bottom part of Fig. 3 shows the charge susceptibility ob¬ 
tained for different values of V together with cxXg. Sur¬ 
prisingly, the agreement is very good even for the largest 
value of E = 0.3, which is already comparable to the 
half-bandwidth D = 1 used in the simulations. 

5. Dynamic Properties 

5 .1 Green’s function 

Top part of Fig. 4 shows the single-particle excitation 
spectrum, —ImG/(£ -I- iS), with M = 5 for several val- 
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Fig. 4. (Color online) Top; Energy dependence of the imaginary 
part of the /-electron Green’s function for channel numbers M = 
5 {main panel) and M = 1 (inset). Bottom: Fit of the spectral 
function by the sum of two Lorentzians at Uf^ = 0.6 for M = 5. 
The parameter values are chosen as V = 0.2, /3 = 200. 


Fig. 5. (Color online) Fitting result for the parameters Ai, Ai 
(top) and A 2 , A 2 (bottom) as a function of with increasing 
channel number M. The parameter values are chosen as F = 0.2, 
/3 = 200 . 


ues of Ufc- As Ufc is increased, a distinct deviation is 
observed compared with the non-interacting lineshape, 
i.e., the Lorentzian; the spectrum exhibits a high-energy 
tail in addition to a sharp peak around £ = 0. Since the 
high-energy tail is not observed in M = 1 as shown in the 
inset of Fig. 4, it is due to the extra screening channels, 
£ 0, without hybridization. 

By analysis of the numerical data, we found that the 
peculiar spectra observed for M > 1 can be well approx¬ 
imated by a sum of two Lorentzians: 


G{appro.)^z) 


Al 

z -b iAi 


+ 


A2 

2: -b iA2 


(16) 


Here, the first term describes the / level peaked at 
£/ = 0, while the second term gives an account of the 
high-energy tail. A fitting result with Ufc = 0.6 for 
M = 5 is shown in the bottom part of Fig. 4. An ex¬ 
cellent agreement is confirmed. The function in Eq. (16) 
is reduced to a single Lorentzian when A 2 = 0, which 
should take place for (/yj, = 0 or M = 1. 

Numerical results for the fitting parameters Ai, A 2 , 
Ai and A 2 are plotted in Fig. 5 as a function of Ufc 
for several values of M. In these calculations, we have 
worked in the Matsubara frequency domain to avoid in¬ 
accuracy caused by analytic continuations. For M = 1, 
we obtain a single Lorentzian, Ai ps 1 and A 2 = 0, as 
expected. The improper result, Ai > 1, is due to influ¬ 
ence of the finite bandwidth used in the simulations. For 


M > 1, the weight Ai is transfered to A 2 as Ufc in¬ 
creases. The spectrum is finally dominated by the high- 
energy tail (Ai < 0.1) in the region Ufc ^ 0.7 for M = 5. 

The energy scale Ai is equal to Aq = 0.063 at Ufc = 0, 
and increases with increasing Ufc in the weak-coupling 
regime. In the strong-coupling regime, on the other hand, 
Ai turns to decrease for M > 1. In particular, an ex¬ 
trapolation of the data for M = 5 to larger Ufc suggests 
an existence of a critical point characterized by Ai = 0 
around Uf^, ~ 0.9. However, we could not reach the crit¬ 
ical point because of a numerical difficulty in the strong¬ 
coupling regime. The difficulty arises not only from the 
increased computational time with increasing M, but 
also from a decrease of the acceptance rate of the Monte 
Carlo updates, which is similar to the case of the ordi¬ 
nary Anderson model with strong Coulomb repulsion. 

The non-monotonic behavior of Ai is a consequence 
of the competition between the exciton effect coming 
from the hybridizing channel and the Anderson orthogo¬ 
nality effect from the screening channels as discussed in 
Section 2. Since there are no additional screening chan¬ 
nels in the case of M = 1, the low energy scale Ai in¬ 
creases monotonically against Ufc- In contrast to Ai, 
the width A 2 of the additional peak (energy scale of the 
high-energy tail) does not show distinct M dependence 
at least for Ufc < 0.5. 
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Fig. 6. (Color online) Energy dependence of the imaginary part 
of the t-matrix for channel number M = 5 for £ = 0 (top) and 
£ 0 (bottom). Dashed line at zero energy in the left panel 

indicates Friedel’s sum rule as l/(7rpo). The parameter values 
are chosen as F = 0.2, /3 = 200. 




Fig. 7. (Color online) Top; Electrical resistivity for channel num¬ 
ber A7 = 5 as a function of temperature for different values of 
[ffc- Bottom: Low energy scales Ai, and A^ together with 
the scaling result A*^. given by Eq. (7) for M = 1 (main panel) 
and M = 5 (inset). 


5.2 t-matrix and transport 

The t-matrix includes all informations about the scat¬ 
tering of the conduction electrons by the local charge. 
Figure 6 shows the C//c-dependence of —Iratii^e + iS) with 
M = 5. We find that Imtf(e) for the scattering channel, 
£ 0, vanishes at the Fermi level {e = ep = 0), while that 

for the hybridizing channel, £ = 0, satisfies the Friedel’s 
sum rule, Imti^o^e = 0) = — l/(7rpo) for arbitrary values 
of Ufc and M. This means that the phase shift S for the 
resonant scattering is tt/2. 

The conduction electron density of states is expressed 
as 

PcA^) = 

TT 

=--lm[g(^z) + giz)tii^z)g{z)] (17) 

TT 

through the £-matrix. Using the properties £^=o(0) = 
—i/{npo) and ti^Q^O) = 0 together with the approxi¬ 
mation g{0) = —inpo, we find that Pc,f=o(0) = 0, i.e. 
the local density of states of conduction electrons van¬ 
ishes at the Fermi level for the hybridizing channel, while 
Pc,r/o(0) = Po: be. pc is unchanged for the screening 
channels. 

To obtain more information about the low-energy dy¬ 
namics, we calculate the electrical resistivity R{T) at ar¬ 


bitrary temperature as 



(18) 


in the Boltzmann equation approach, where the relax¬ 
ation time T(e) is obtained from the t-matrix in the hy¬ 
bridizing channel^^^ as T(e)“^ = —2lm.ti=o{e -I- iS). The 
temperature dependence of R{T) for M = 5 is shown in 
the top part of Fig. 7 for several values of Ufc. Here, the 
temperature is scaled by a characteristic energy Ap (see 
below for detail). As an accuracy check of the numerical 
data, we also show the analytically derived resistivity for 
Ufc = 0 with ti^Q{z) = V'^Gf{z). We confirm from the 
numerical data that the scaled resistivity R{T) / R{0) is 
a universal function of T/Ap in the temperature range 
shown in this figure. 

At low temperatures the resistivity behaves as 


Rlfl{T) = Ro 



(19) 


within the local Fermi-liquid theory. Using this expres¬ 
sion, we have determined the energy scale Ap. Here, we 
assumed that Rq and a are independent oiUfc, which can 
be determined from the Sommerfeld expansion of R{T) 
as Rq = 2Q:/(7rpo) and a = air^jd. The result obtained 
for A/j is shown in the bottom part of Fig. 7 for M = 1 
and M = 5 together with the low energy scales Ai and 
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Fig. 8. (Color online) Top: Energy dependence of the numerically 
obtained /-electron self-energy and parameter A for M = 5 at 
Ufc = 0.6. The parameter values are chosen as V = 0.2, (3 = 200. 
Dashed lines in the low-energy range correspond to the quasi¬ 
particle results to be given later in Eqs. (33) and (34), while 
dashed green lines show the approximation given in Eq. (21). 
Bottom: Wave-function renormalization factor af for M = 5 as 
a function of Ufc- Dashed green line corresponds to the quasi¬ 
particle perturbation result ay = 1 — (u/uc)^ to be given in 
Eq. (35). 


= 27rT[tp~^(27r^Txoi'r)) — 1/2] of the local Green’s 
function and charge susceptibility, respectively. We find 
that the energy scale A/j coincides with Ai. Thus we con¬ 
clude that the low-energy scale for the t-matrix matches 
with the low-energy scale of the local Green’s function. 

5.3 Self-energy 

We now discuss the self-energy to get more information 
on the peculiar spectra of Gf{e). Top part of Fig. 8 shows 
the energy dependence of the numerically obtained local 
self-energy S/(e) together with A(e) = Trpo\V -|-S/c(e)P 
for M = 5 under finite Ufc- While S/ shows substantial 
energy dependence, the quantity A involving the self¬ 
energy E fc is almost independent of energy. A crucial dif¬ 
ference with the ordinary Fermi-liquid behavior is shown 
in the high-energy region, ImE/(£) ^ const. It is related 
to the high-energy tail of ImG/(e). 

We find that the numerical data for ImE/(e) is well 
approximated by the formula 

= ( 20 ) 

where we use the convention that F, Ag > 0. The 
Kramers-Krdnig relation constrains the complex self¬ 


energy to be 


E^"PP’'°x)(£) 


eF 


e -I- iAs 


( 21 ) 


Fit with the approximation (21) for both the imaginary 
and real parts of Sy is shown in top part of Fig. 8 by 
dashed green line. Formula (21) reproduces the Fermi- 
liquid properties ImEy(e) oc —and ReSy(e) oc —e in 
the low-energy range jej < (!l(As). We note that the be¬ 
havior ImE/(e) ~ const, at large energies should break 
down for £ ^ The apparent constant behavior sug¬ 
gests the presence of an additional characteristic energy 
much larger than Ag. 

The U /c-dependence of the wave-function renormaliza¬ 
tion factor a/ for M = 5 is shown in the bottom part of 
Fig. 8, which we obtain from the numerical self-energy 
data as 


a/= [1 - ImS/(i£o)/£o] (22) 


The decreasing feature of a/ with M > 1 as the Goulomb 
interaction is increased indicates the formation of a corre¬ 
lated Fermi-liquid state with large quasi-particle effective 
mass m* since m* ~ l/o/- Vanishing of a/ indicates the 
quantum critical point where the local charge decouples 
from the conduction electrons. Unfortunately, this criti¬ 
cal point is difficult to be reached by the continuous-time 
quantum Monte Garlo method because of the numerical 
difficulties mentioned in Section 5.1. 


5.f Relation between energy scales 

So far, we derived two sets of energy scales from the 
numerical data: {Ai, A 2 } from the Green’s function, and 
{F, Ag} from the self-energy. Those are related with each 
other. In the following, we derive the formula connecting 
them. 

From Eq. (A-2), the /-electron Green’s function G/(£) 
is given by 

^ £ + *A-E/(£)' 

Here, A denotes an effective hybridization strength de¬ 
fined by A = 7rpo|U + S/c(0)p. We neglected T,c,e=o and 
the energy dependence of E/c(e) according to the nu¬ 
merical results shown in Fig. 8. By replacing T,f with 
Ej^PProx) given by Eq. (21), and equating G/ to the two- 

Lorentzian form qGpp'^°^'> given by Eq. (16), we obtain 
the following relations: 


A 1 A 2 -|- (1 — Ai)Ai — Ag, 

A 1 A 2 = AAg, 

Ai-b A2 = A-f Ag-b F, 
A 2 = 1 — Ai. 

Solving Eqs. (24)-(27) for Ai and A 2 we obtain 
Ai^2 = 2 ^ ^'i) 

T v^(r + A + Ag)2 - 4AAg' . 


(24) 

(25) 

(26) 
(27) 


(28) 


Now we compare {Ai,A 2 } with other energy scales 
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Fig. 9. (Color online) Energy scales Ai (top) and A 2 {bottom) 
together with the energy scales As and A as a function of Ufc for 
M = 5. The parameter values are chosen as y = 0.2, (3 = 200. 
Dashed green lines correspond to the quasi-particle perturbation 
result to be given later in Eq. (42). 


such as A. Figure 9 shows their Ufc dependence for M = 
5. The enhancement of A over the bare hybridization is 
due to the exciton effect. Concerning Ai 2 , we find that 

(i) Ai shows a crossover from A to As/2 with increasing 
values of Ufc, while 

(ii) A 2 interpolates between A and 2A. 

We will discuss these characteristics by a microscopic 
Fermi-liquid theory in the next section. 


6. Quasi-Particle Perturbation Theory 

6.1 Second-order self-energy 
For descriptions of low-energy properties, we may work 
with a quasi-particle Green’s function. Assuming the or¬ 
dinary Fermi-liquid properties for T,f{e) in Eq. (23), we 
obtain 


G}{e) 


1 

e + iA*-S}(e)’ 


(29) 


where Gf = afGf A* = a/A, and Sj = af'Sf with 


Of = [1-dT,f{e)/de]^^g, (30) 


being the wave-function renormalization factor. Since 
Of < 1, G*j, has a smaller characteristic energy scale, 
A* < A, and quasi-particle descriptions are valid in the 
region jej < 0{A*). 

We now derive S)- by a perturbation theory with re¬ 
spect to Ufc- The leading contribution is given by the 
second-order diagram shown in Fig. 10. The explicit ex¬ 


es 



Fig. 10. (Color online) Second-order self-energy diagram with 
dashed line corresponding to local (/) state and solid line repre¬ 
senting conduction electron state. Dots represent the Coulomb 
interaction Ufc- 


pression reads 

711712 ^ 

where Qc,e. is the local component of conduction-electron 
Green’s function renormalized by the hybridization pro¬ 
cess. The important point here is that the influence of 
Ufc in the hybridizing channel £ = 0 is already included 
into G*j as the effective hybridization A. Since the lo¬ 
cal density of states Pc,^=o(0) vanishes at the the Fermi 
level as presented in Sec. 5.2, the channels contributing 
to —ImS^(£) are those with I Q. The Green’s functions 
Gc4=jto of these non-hybridizing channels may be replaced 
by 5fc(fe) = (1/V) ~ £k)~^■ Performing the Mat- 

subara summations and taking the imaginary part, we 
obtain 

nOO pQ nOO 

-lmE}{e) = u^{M-l) dei de 2 dss 
Jo J — 00 Jo 

x27rpj(ei)(5(e - £ 1 - 1-62 - £ 3 ), (32) 

where p). = —(l/7r)ImG)'. In the low-energy limit, £ <C 
A*, we may replace G*j:{ie) by G*jP(ie) = (ie — iA*)“^, 
and the integrals can be performed analytically to yield 

=u^{M-l)e‘^/A* (33) 

with neglect of terms of order A*(£/A*)^. 

The real part of S)- may be evaluated using Eq. (31) as 
well. However, in contrast to the imaginary part, high- 
energy processes give finite contributions in this case. 
Then we cannot use the low-energy form for Gf We in¬ 
stead derive the low-energy limit of ReS). by the follow¬ 
ing argument. The particle-hole symmetry ensures the 
condition ReE)-(O) = 0. It follows that the low-energy 
expression for ReEy may be given within 0{v?) as 

ReE/(£) = -{ulucfe, (34) 

where Uc is the critical value of poUfc that gives the 
quantum critical point. According to the second-order 
scaling we obtain Eq. (13) which reduces to —>■ 1/M 
in the limit of M ^ 1. Since <C 1 for large M, the 
second-order self-energy should give precise account for 
u < Uc- We then extrapolate the second-order theory 
for smaller M by modifying the form to = 1/(M — 
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1), which avoids correctly the critical point at M = 1. 
Combining with Eq.(30), we then obtain 

a/ = 1 — (M — 1 )m^ = 1 — {u/ucf ', (35) 

which is reduced to zero at it = Uc = {M — Fi¬ 

nally, the quasi-particle Green’s function is obtained from 
Eq. (29) as 

G){e) = {e + + i{l - af)e^//S*Y^ (36) 

by using Eqs. (33), (34), and (35). 


Namely, we obtain from Eq.(28) the following relation: 


A 1 . 2 /A = 1 =F ^1 - a/ = 1 =F u/uc (42) 

by using Eqs. (35), (39), and (40). The limiting behavior 
is thus given by 


Ai/A = l- 


and 


a//2 = A,/(2A) ifa/<l 
1 if a/- Ri 1 ’ 

(43) 


6.2 Comparison with numerical data 

Comparing the approximate formula (21) of the self¬ 
energy with Eqs. (33) and (34), we obtain the following 
correspondences between the phenomenology and quasi¬ 
particle theory: 


5ReE/(e) 


de 


£ = 0 


1 i9^ImEj(e) 


2 e=o 

which give the relations 


T/As = aY - 1 , 

r/Al = {aY-l)/A*, 


(37) 

(38) 


A, = A*=afA, (39) 

T = {l-af)A = {M-l)u^A. (40) 

Thus, the phenomenological parameters T and A^ in the 
approximate self-energy are determined by the single pa¬ 
rameter Of for given A. This is a strong constraint im¬ 
posed by the quasi-particle Fermi-liquid theory. It is thus 
interesting to check the accuracy of the quasi-particle 
perturbation theory in the light of the numerical results. 

First we examine the relation given in Eq. (35) which 
expresses the wave-function renormalization factor af in 
a simple way. The G/c-dependence of the estimate 1 — 
{u/ucf = 1 —(M—l)u^ with M = poUfc is shown as green 
dashed line in the bottom part of Fig. 8 together with 
the numerically obtained aj. We can observe that the 
simple, second-order quasi-particle formula describes the 
wave-function renormalization factor reasonably well. 

Next we discuss the case of self-energy. Top part of 
Fig. 8 includes also the fit of the numerical self-energy by 
the result of the quasi-particle theory given in Eqs. (33) 
and (34) as dashed lines. We conclude that the theory 
works well in the Fermi-liquid range |e| < A*, i.e. at low 
energies. 

Using the relations (39) and (40), the approximate self¬ 
energy given in Eq. (21) is expressed as 


^(approx), X _ ^ (1 af)Ae 

^ e + iafA 


(41) 


in the quasi-particle theory. This formula is also shown 
in the top part of Fig. 8 by green lines, and gives an 
excellent fit of the numerical self-energy in the whole 
energy range. This is not surprising since the approx¬ 
imation given in Eq. (21) has only two parameters, A^ 
and r, and if the quasi-particle theory fits the self-energy 
around e ~ 0, it will fit also the curve in the whole energy 
range. 

Finally, we discuss the energy scales Ai and A 2 . 


A2/A = 1 + vT^= (44) 

(1 if a/ « 1 

Actually, these are exactly the limits that we found from 
the analysis of the numerical data (see Fig. 9). We show 
the expressions A(1 =F u/uc) from Eq. (42) in Fig. 9 by 
dashed green curves as well, which describe Ai and A 2 
relatively well. We plotted the energy scales Ai and A 2 
in Fig. 9 instead of the ratios Ai ^/A because numerical 
errors are larger in Ai^ 2 /A in the low-{7/c range since the 
fit with two Lorentzians is unambiguous in this region; 
the contribution from the wider Lorentzian is very small. 

We conclude that the extrapolated quasi-particle the¬ 
ory describes the U/c dependence of both the energy 
scales {Ai, A 2 } and the self-energy in a surprisingly wide 
range of Ufc. However, this theory does not give any 
microscopic origin neither for the constant behavior in 
ImEy(£) at large energies nor for the wider Lorentzian 
component in Im G/(e) since these properties appear out 
of the valid energy range of the perturbation theory. 

7. Summary 

In this paper we have studied the multichannel inter¬ 
acting resonant level model by means of the continuous¬ 
time quantum Monte Carlo method in a wide range of 
the Coulomb interaction U/c and channel number M. 
Thermodynamic and dynamic properties have been de¬ 
rived accurately and have been discussed in the light of 
analytic approaches. 

We find that thermodynamics of MIRLM such as the 
local charge susceptibility is entirely determined by the 
single energy scale A*^ within the scaling theory. On the 
other hand, dynamics contains multiple energy scales be¬ 
yond the simple scaling theory. For example, we find that 
the single-particle excitation spectrum acquires a high- 
energy tail for M > 1 with increasing values of Ufc in 
addition to the narrow resonance around the Fermi level. 
This composite line shape can be well described by the 
sum of two Lorentzians, where the narrow one with en¬ 
ergy scale Ai corresponds to the / level peaked at e/ = 0, 
while the wider one with the scale A 2 gives an account 
of the high-energy tail. We find that while the narrow 
energy scale Ai shows strong M-dependence and non¬ 
monotonic behavior for M > 1, the larger energy scale 
A 2 is independent of M for moderate values of Ufc. 

The numerically obtained self-energy also shows un¬ 
usual behavior: its imaginary part is nearly constant in 
the high-energy region, which is related to the tail of the 
Green’s function. We note that a three-peak structure 
similar to the case of the symmetric Anderson model is 
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expected for Im E f close to the quantum critical point 
since the renormalized hybridization tends to zero, and 
therefore addition or removal of an / electron from the 
ground state may accompany the additional peaks. Ac¬ 
tually, preliminary numerical data showing this situation 
has already been obtained numerically.^^^ 

A quasi-particle perturbation theory from the Fermi- 
liquid fixed point is used for the microscopic understand¬ 
ing and description of the low-energy part of the single¬ 
particle spectra. The microscopic theory provides a con¬ 
straint among the parameters in the phenomenological 
theory, and gives description of the spectrum by a single 
parameter a/ for given A. 

Finally, we propose that the multichannel interacting 
resonant level scenario might be responsible for the pecu¬ 
liar heavy-fermion state of certain Samarium compounds 
with large mass enhancement and magnetic field insensi¬ 
tivity. Namely, we speculate that in the regime near the 
quantum critical points u± with vanishingly small effec¬ 
tive hybridization a highly renormalized heavy-fermion 
state is developed with large effective mass, which is inde¬ 
pendent of the external magnetic field since only charge 
is involved. It remains to see how large is the number of 
active channels in real materials. 
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Appendix A: Renormalized Green’s functions 
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Viz)Gf{z) 
g{z)-^ - Ec,o(^) 

GcA^) = g{z)+g{z) 
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g{z)Y.cAA 
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